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Abstract 

This paper is devoted to the use of the entropy and duahty methods for the existence theory 
of reaction-cross diffusion systems consisting of two equations, in any dimension of space. Those 
systems appear in population dynamics when the diffusion rates of individuals of two species depend 
on the concentration of individuals of the same species (self-diffusion), or of the other species (cross 
diffusion). 

1 Introduction 

We are interested in populations consisting of individuals belonging to two distinct species (or two classes in the same 
species) of (typically) animals, and interacting through diffusion and competition. 

The unknowns are the concentration (number density) of individuals of the first species ni := ui{t,x) > and of the 
second species U2 := U2{t,x) > 0. 

In absence of competition, the respective populations ui and U2 would grow at a (strictly positive) rate ri and r2. 

The competition is taken into account through logistic-type terms, in such a way that the growth rate becomes 
ri — sii(ni) — 512(112) for the first species, and r2 — S22(m2) — S2i(«i) for the second one. The (nonnegative) terms sii(ni) 
and 522(1^2) arc called intraspecific competition, while the (nonnegative) terms 812(^2) and S2l(wi) are by definition the 
interspecific competition. Since we are interested here in situations when there is no cooperation (or predator-prey type 
interaction), we shall assume in the sequel that all functions Sij are nonnegative and increasing (see 1221 . p. 94). 

We also assume that the individuals of the two species diffuse with a rate that depends on the number density of both 
species: we denote dii{ui) + ai2{u2) the diffusion rate of the first species, and d22{u2) + <J2i(mi) the difi^usion rate of 
the second species. The difi^usion terms related to (iii(ui) and ^22 ("2) are called self diffusion, those related to ai2(ji2) 
and 021 (mi) are called cross diffusion. We are interested in this paper in the case when all the functions da, aij are 
nondecreasing (that is, all individuals try to leave the points where the competition is highest). This type of model were 
introduced in 1241 . 

In the sequel, we systematically denote 

a.ii{ui) = Ui dii(ui). (1) 

We are led to write down the following system (on R_|. X f2, where SI is a smooth bounded open subset of with outward 
unit normal vector n{x) at point x S dO.): 

9t«i - A[aii(«i) -I- uiai2(M2)] = fi (ri - sii(mi) - Si2(n2)) := -Ri2('!ii, M2), (2) 

dtU2 - A [022 ("2) + «2a2l(wi)] = 1*2 ('•2 - S22(W2) - S2l(«l)) := -R2l(w2,'"l), (3) 

with the homogeneous Neumann boundary conditions 

Vmi ■ n{x) = Vu2 ■ n{x) =0 on R+ X dVl, (4) 



and the initial condition 

ni(0,-) ="1" > 0, W2(0, ■)=«!"> 0. 



(5) 



The case treated in this paper corresponds to a situation where the reaction terms are strictly subquadratic or dominated 
by the self diffusion, and where the cross diffusions are subquadratic. Extensions to other cases with faster growth of 
reaction terms or cross diffusion pressure will be studied in a future work. 

We detail below the set of mathematical assumptions that will be imposed on the coefficients a^j , Sij, and which 
correspond to the case described above: 

HI For i,j £ {f,2}, S IR+ and Sij is a nonnegativc continuous function on R+ which is either strictly sublinear: 

lima:_>+oo "^^^^^ = 0; or dominated by the self diffusion in the following sense; lim^j—j+oo xd^^l) ~ ^ ^^'^ (^'-'^ * ^ ■5') 

H2 For i j £ {1,2}, a^j is continuous on I84- and belongs to V^QO, +oo[). It is also nonnegativc, nondccrcasing, 
concave and vanishes at point 0. Furthermore, there exists o G]0, 1[, C > 0, such that 

Vx 6]0, +oo[, x^a'^jix) > C. 

H3 The self diffusion rate da (recall notation JTJ) is continuous on R_|_ and belongs to '^^-'^{JO, +c«[). It is also nonnegative, 
nondecreasing and such that dii(O) > (note in particular that a'^^ is bounded below by a strictly positive constant). 
In the proof of our main Theorem, we use da > 1 for the sake of readability. 

The union of the assumptions HI, H2 and H3 on the parameters [for i,j S {1, 2}] r^, Sij, aij (and thus da) will be called 
in the sequel the H assumptions. 

Since we consider the homogeneous Neumann boundary condition, it is useful to introduce the following notation : 

Definition 1.1. for any space of functions defined on Q whose gradient has a well-defined trace on the boundary dQ 
(such as '(f°°(r2) or ll^{Q) for instance), we add the subscript v (the former spaces hence become '(f„°°(r2) or Hj(r2)J to 
describe the subspace of functions satisfying the homogeneous Neumann boundary condition. 



Our main Theorem reads: 



Theorem 1.2. Let ft be a smooth (6"^) bounded open subset ofR'^ (N >1). Let := {u^^,ui^) e L^(Q.)'^ be a couple 
of nonnegative functions and assume that assumptions H are satisfied on the coefficients of the system. 

Then, for any T > 0, there exists a couple u := {ui,U2) of nonnegative weak solutions to ^ - Q on [0, T] in 

the following sense: for i = 1,2 (j i), 



ri 

Jo Jn 



dii(ui{t,x)) + aij{uj{t,x)) 



\ui{t, x)\'^ dx dt < +00, 



)dx dt 



and for any 9 S T\;,V()° (Q,)) , we have the weak formulation 

n roc n 

- / u^{'(x)e{fi,x)dx- / / u^{t,x)dt6{t,x) 
Jn Jo Jn 

-J A6'(t, x) l^aii (mi [t, x)) + ui (t, x) ai2 (u2 {t, 

= / / Ri2{ui{t, x),U2{t, x)) dx dt, 
Jo Jn 



dx dt 



and 



) dx dt 



- f ui^'-(x)9{0,x)dx - f f U2{t,x)dte{t,x)< 
Jn Jo Jn 

^ j^^ j A6'(t, x) |^a22('ii2(t, 2^)) + «2(i, a;)a2i(ni(t,x))j dxdt 
= / / R2l{u\{t,x),U2{t,x))i. 

Jo Jn 



) dx dt. 



Moreover, denoting Qt ■= [0, T] X Q, the following bounds hold: 



1 



a2i('Ui)V'Ui 



ui 

Vv'a2l(ni) 



U2 



-a'j^2('^2)V«2 



<^T(l + ||u«||i. ), 



L^(Qt) 



+ Vv'«12(w2) 



+ Vv'a2i(«i)ai2(w2) 



for some positive constant Kx depending only on T and the data of the eguations (oij, ri, Sij). 



(6) 



< Kril + \\u° 



^L\n)> 



Note that all the terms in the weak formulation above are well-defined thanks to assumptions H. Indeed, remembering 
that d-ii is bounded below, estimate ^ implies Ui 6 L'^([0, T] X Q) so that the reaction terms arc intcgrablc (thanks to 
HI). Using the identity 



S12(«2) ui 



812(1*2) 



\/d22iui) + ai2{u2) 



«2 v''^22(«2) + ai2(n2) \ 



^012(^2) 



- 



ui v' ai2(n2), 



the terms coming out of cross diffusion are also well defined due to concavity assumptions H2. Those coming out of self 
diffusion are intcgrablc thanks to 
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Wc now comment the assumptions on the initial data and coefficients of system ((2]l-|[3}. The requirement that S 
(L^(f!))^ can certainly be relaxed to it" g H~^(S7)^ and tpi{u'-") S L^(f2), with tpi defined later in Definition 12.11 Cthis 
should be enough to keep estimate ||6ll). It is most likely possible to relax assumption HI, provided that H2 and H3 are 
reinforced. Finally, there is some hope of treating the special case when reaction terms are exactly quadratic, thanks to 
recent improvements in the theory of duality Lemmas [5]. 

We restricted in this paper our study to the case when both a^j , i ^ j are concave, whereas our methods should adapt 
in the case when one of them is concave and the other one is convex. The corresponding theory is then quite different 
and will be left to a future work. Note that when both a^j, i ^ j are convex, our feeling is that existence of weak global 
solutions does not hold in general. 

Our opinion is that for systems involving cross diffusion and consisting of more than two equations, the type of 
Lyapounov functional that wc build in the sequel can exist only for a very small class of cross diffusion terms (that is, 
strong algebraic constraints have to be assumed on the cross diffusion coefficients). 

Finally, we do not treat the case when one cross diffusion term is missing (sometimes called the "triangular case"), 
since the Lyapounov functional that we introduce in the sequel degenerates in this case. 

Let us also explain the meaning of assumptions H when all the functions appearing in Q-JSj arc of the form x i— > 
(or X I— > X + x' for the self diffusion), with (j E M. In this case, if Sij(x) = Sij x'^'J and aijix) = Di x&ij + Aij x">J , with 
Sij , Aij , Di > and crij,aij £ R, those assumptions become: 

V't=l,2, < o-ii < sup(l, Oii); 

ViT^i, < < sup (^ "''^ \ 1 + ^ ^ , 
< aij < 1. 

As a consequence, our Theorem provides existence of global weak solutions for systems like 

dtui - A[(Di + All ui + Ai2 M'2 ) ui] = «i(ri - 5ii «i - S12 "2), (7) 
dtU2 - A[{D2 + A21 uj + A22 "2) "2] = "2 ('"2 - •921 ui - S22 U2), (8) 
with < «,(5 < 1; /3,7 > 0, and Di, Aij , Sij > 0. 

Let us now describe how our work fits in the existing literature. 

The question of local and global existence of classical solutions of systems like (l2]l-|[3j was treated in many particular 
cases. Most of them deal either with the case when the system is in fact parabolic (that is, when the diffusion matrix 
is elliptic), which amounts to assume that self diffusion is dominant w.r.t. cross diffusion (for instance when 8 An > 
j4i2,8j422 > A21 in eq. Q, ^ with o = /3 = 7 = <5 = 1, cf. [26]), or in the triangular case (when a2i(ni) = 0, cf. 1201 '). 
The question of local existence is usually treated using Amann's theorem [l], and extension to global existence requires 
additional structure to preserve boundedness of solutions see [8] 1121 for instance. 

Our work is concerned with situations that are neither parabolic, nor triangular. It extends to a very large class 
of problems the result of [7], dedicated to the cross diffusion model for population dynamics introduced by Shigesada, 
Kawasaki and Teramoto in I24| . where the cross diffusions have a linear form, that is /3 = 7 = 1 in eq. (0, (|8]l. Note that 
this system can be seen as a limiting case of the equations that we treat. 

Our results rely on two main ingredients: entropy structure and duality Lemmas. We show that our systems possess 
a hidden entropy-like structure, strongly reminiscent of the entropy structure exhibited in [7]. In general, this structure 
however gives less estimates than in [7]. We therefore need another ingredient in order to recover existence of weak global 
solutions, namely duality Lemmas: we recall that duality Lemmas enable to recover L^ type estimate for solutions to 
linear singular parabolic equations (with variable coefficients) when the diffusion rate is inside the Laplacian. This is how 
estimate (|6]l is derived. 

For the use of an underlying entropy structure, its link with symmetrization, and its applications to existence of weak 
solutions, we refer to 1181 [9l . The possibility to use such a structure in the case of cross diffusion was first noticed in 1131 . 
and exploited in [7] [G] 1171 1151 . The duality estimate that we use comes from 1231 , and was applied together with entropy 
methods in the framework of reaction-diffusion systems in 1111 . Is was also applied in the framework of cross diffusion or 
similar models in [2] El HI- 

We finally quote some works dealing with other aspects of cross diffusion models. For modeling issues, we refer to 
|16l 1101 [3I [2I. For the analysis of equilibria, we cite 1161 1191 for instance. 

Unfortunately, as often in papers dealing with cross diffusion, the process of approximation enabling to make the 
estimates and structures rigorous is quite involved, and gives rise to various difficulties which explain the length and 
technicality of the proofs. 

Our paper is structured as follows: in Section[2] we introduce notations which are used in the proof of our main Theorem, 
especially those related to the Lyapounov functional that we systematically use. We also present some classical lemmas 
used in the sequel. Then, Section \3\ is devoted to the proof of existence of a solution to a finite-dimensional (Galerkin) 
approximation of a discrete time version of our (smoothed) system. A priori estimates (and their dependence with respect 
to the various approximations) are provided for such solutions. We let the dimension of the Galerkin approximation go to 
infinity in Section U] The duality estimate is presented and proven in Section [5] The last section (that is, Section [5} is 
devoted to the relaxation of all remaining approximations, which leads to the proof of our main Theorem. 
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2 Preliminaries 



2.1 Entropy structure 

We first introduce some notations which enable to rewrite our system under a form in which a Lyapounov functional 
naturally appears. 



Definition 2.1. For given cross diffusion parameters ai2 and 021 (satisfying assumption U2), we introduce ipi, ipi (for 
i = 1,2), as: 

a'(t) , a'(t) 



tpi(x) := <fiiit)dt. 



We also define 

Wi{t, x) := ipi{ui{t, x)). 
One can then rewrite the system in a symmetric form: 

~A{ui,U2) 



/ a-i -I (uT )-\-ai2(^2) \ 

The terms Vui^ have to be considered as scalars for the matrix product, the divergence being understood line by line (after 
the matrix product). 



Multiplying the equations of the system by wi , U12 and integrating in space we obtain formally the entropy identity 

d 
dt 



■ JjM^l)+i>2{u2)) + Jj\'wi,\'W2)A{uuU2) (^Zl) = ("'l''"2) (^^21) ■ 



2.2 Properties of the symmetric matrix 

For ui,U2 > 0, u := (mi,«2), the matrix A{u) is defined by 

A{ui,U2) := "^i*"!^ , , , , , (11) 

and we denote the associate quadratic form Q(u). For a function u = {ui,U2) : H — >]0,+oo[^, we defined in the previous 
subsection w = {w\,W2) = ('Pi('i'i), ¥'2("2))- We will have to deal with the expression 

Q{u){Ww), 

since this term naturally appears in the entropy estimate. We establish the following proposition, that will be useful later. 



Proposition 2.2. Under the H assumptions, the application A : ]0, +oop — > A^2(IR) (defined by (j j|) ) belongs to 
■^"(JO, +oop) and takes its values in S^^(R) (space of strictly positive symmetric matrices). We have furthermore the 
two following estimates on the quadratic form Q 

Q{u){^w) > — a;,i(Mi)|Vni|2 + — a'l2(«2)|V^t2|^ (12) 
111 U2 

Q(u)(Vii)) > 4| |vVa2i(Mi)|^ + |vVai2(«2)|^ + | V Va2i(Mi)ai2(«2)p }. (13) 



Proof. The fact that A £ '^''(JO, +oo[2) is an easy consequence of the regularity of the functions aij. As for the strict 
positiveness of the matrix, we just decompose A between self (clearly strictly positive) and cross diffusion: 



^(ni,«2)=r^ ' + 



4^«2/ I «1«2 f?4S)"2 

£il2("2) / \ fj2lU2j 

We see that due to the assumptions on ai2 and 021, the cross diffusion matrix C'(ui,U2) is nonnegative. Indeed, for 
U1,U2 > 

0-12 ('"2) 121 («l) . 

i~/ — r"i ^ ~r~/ — 7'"2 > "i"2 X U1M2 

a'.2i(ui) «12("2) 

t 

0,12(^2)0,21 (ui) > uiU2a'i2(u2)a2l("l)i 
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which is true by concavity. Wc get 



Q(u){Vw) = [*Vii;] [A{m)Vio] > [*Vii;] [B{u)\'w] = °"|"^| mi|Vwi|^ + °22("2) ^|y^^|2^ 



Since Wi = ipi{ui), with </5^(a;) = — , we end up with 

X 

Ul U2 

which, using assumption H3, leads to the first lower bound I I12II . On the other hand, expanding directly Q{u){Ww) with 
the definitions of A{u) and w, we get 

Q{u){Vu,) = KiK) + ai2(«2)] ^2i(!^|v«i|2 + [a'^^{u2) + a2i (m)] I V«2 P 

+ 2 a2i(-ui)ai2("2)(V'Ui, Vm2)- 
The concavity property aji(ui) > a'^^{ui)Ui used before can be written (here for i = 1, j = 2) 

"21 ("i) > "21 ("i)^ 
Ml ~ (121 (mi) 

and together wc the assumption H3, we get 

Q{u)(V«,)>aii(«i)^^^|V«ip + 42{n2)^i4^l^"2|^ 
a2i(«i) ai2(n2) 

+ ai2(»2) "^'^'"\ |Vtti|^ + a2i(m) "'^^"^\ |Vn2|^ + 24i(Mi)a'i2(M2)(Vm,VM2) 
021 (Ml) ai2('li2) 

= 4| VvWKTP + 4|VV«12(«2)|^ 



+ 4|Vai2(«2)VVa2l(Ml)P + 4|i/a2l(«i)VVai2(«2)|^ 



+ 4 X 2v'a2i(«i)v''ii2(«2)(Vv'a2i("i), V 1/112 (m2)>, 
and we thus end up with eq. II13II . □ 

2.3 Properties of the functions (fi and 

We shall need in the sequel the following elementary result: 

Lemma 2.3. Take h,l S 'i^'''(R+) fl '^-'-(JO, +oo[) with h concave and I nonnegative and convex, with i'{x) > for all x 
large enough. 

Then there exists a constant Af^ £ > such that h(x) < Af^ i (1 + £{x)) for all x 6 R+. 

Proof. If h is bounded from above, then Af^ i = sup/i works. Otherwise, /i' > on R+. Say that < i'{x) for x > M. 
Then for all x&R+,x-M< Then since h > and h is concave, wc can write 

h{x) < h{M) + h'{M)(x - M) < h{M) + - m)) < KM) + Y^^i.^), 

so that the constant Af^ i = max (ji{AI), jtjjjj^ works. That concludes the Proof of Lcmma l2.3l □ 

Using the previous assumptions on the Oij, one can see that the ipi,il>i at least belong to'i#'^(]0, +oo[). The ipi are 
convex functions and the ipi strictly nondecreasing. Moreover: 

Lemma 2.4. We assume H2 on the coefficients Oij (i ^ j). Then the ipi are convex functions and the ipi are strictly 
nondecreasing. Moreover (for i ^ j): 

(i) For all x,y £ R+, i)'^{x){x - y) > '4)i{x) - ipi{y). 

(ii) ^[{x) = o^_jQ+ (1/x) and hence x'^l)'^{x) > B, for some constant B < for all x £]0, +oo[. 
(Hi) ipi has a limit at point 0"*" fV'l(O) = ci2l(l)j, furthermore ipi is strictly positive on R_|_. 
(iv) There exists a constant D > such that, for all x S R+ 

Vae[0,l], x" +aji{x) < D{l + ^Pi{x)), 
x^'^(x) < D(l + iP,[x)). 

Proof. Let us treat only the case i = 1, the other one being similar. 

(i) ipi is convex. 

(ii) tp'ii^) = 'Pli^) a-iid 

, s C"^ a'yJt) ra2i(t)l^ r a'2i(t) , , , 

Jl t L t Ji f2 x^0+ 

since 021 S '^''(R+) and 021 (0) = 0. The function x 1— > Xipi(x) is strictly nondecreasing after x = 1 and bounded 
near 0, hence lower bounded. 
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(iii) Note that V" (*)* = CjiC*)- We have 

a2l{l) - a2i(x) = J #i'(t)dt = [tVi(t)]^ - y i''i{t)dt, 

then noticing that = = by construction, we have 

a2i(l) - a2i{x) = tpi{x) - Xipi{x). 

hence the previous point (ii) gives the limit near 0. For the positiveness, just notice that ipi = ip'- is negative on 
[0,1]. 

(iv) We use Lemma 12. 31 with h{x) = x + aji{x), l(x) = ^i{x) and x" < 1 + x in order to obtain the first inequahty. For 
the second inequahty, we use the same Lemma with h{x) = x^'^^{x) — 2^i{x), and £{x) = tpi{x). 

□ 

2.4 A small perturbation 

Since ipi may not be one to one, we use a small perturbation of this function and consider the following definition: 

Definition 2.5. Let us assume H2 (and recall Definition \2. Id on the coefficients aij (i j^j), and H3 on the coefficients 
an . 

For all £ > (small enough), we introduce 

flix) := ipi{x) + £ln{x), 

and, equivalently , (for i 7^ j) a|^. := a^j + ex, and 4'i{x) '■= ^i{x) + £a:ln{x) — ex. 

We also introduce = with d?^ = -(^{dii), where (7e)£>o is an increasing family of smooth nonnegative 

and nondecreasing functions, such that 7^ < e^^ on R+, 7e(l) = 1 and {7e)e>o is uniformly converging to the identity on 
compact sets. 

Finally, we denote by A'^ the matrix A (defined by ). where the coefficients aij are replaced by (for all 
i,3 e {1,2};. 

The crucial point is the following: Proposition 12.21 is still true when one replaces the functions a^j and A by their 
£-approximations, and A^, (assumption H2 and H3 hold for the afj), and if one tries to reproduce Lemma 12.41 with 
these new functions, all the inequalities remain the same, the constants being a little bit changed but not depending on e. 
We write the following lemma, which summarizes the situation: 

Lemma 2.6. We assume that H2 and H3 holds on the coefficients Oij (for all i,j G {1,2}), and use the notations of 
Definition \2.5\ 

Then, the coefficients a?^. also satisfy H2 and H3 (with constants independent of e for e > small enough). Also B 
and D (the constants of Lemma \2.4\ l may be changed in order to have, for i = 1,2 and < e < 1.' 
(i) For allx,y& R+, (x)(x - y) > (a;) - (s/)- 

(ii) ipf'{x) = o^_jQ+ (l/x) and furthermore, for all x £]0, +oo[, a;')/)|'(x) > B — ee~^ . 

(iii) ipi has a limit at point 0^ . For e small enough, f/)^ is strictly positive on R+. 

(iv) For all :r S K+ 

Va 6 [0, 1], + al^{x) < D{1 + e) (1 + i/>f (x)), 

xV'r(x)<D(l + e){l + ^f(x)), 
where D is the constant defined above in lemma [2.4\ 

Proof. As before, we treat only the case i = 1. 

(i) ipf is still convex. 

(ii) ln(a;) =o^_jq+(1/x) and exln{x) > —ee~^ > —e~^. 

(iii) X ln{x) — x goes to with x and e{x ln{x) — x) > —e. 

(iv) Firstly, we notice that D > l + aji{l) > 1. For the first inequality, is sufficient to notice that since i/if < tpi, we have 

ex < e{x + aij{x)) < eD{l + tpi), 

and we conclude with 

x'^ + alj{x) < D{1 + e){l + i/>,(x)) < D{1 + e)(l + (x). 
For the second inequality, we use 

£a;ln(a;) < e{x\n{x) - x + 1) + ex < De{x\n{x) - x + 1) + D{1 + ipiix)), 

so that 

xtPI'ix) < x^P[{x) + exln{x) < D{1 + ^i{x)) + De{xln{x) - x + 1) + D{1 + ^i{x)) < D{1 + e){l + -iPf {x)) 

for some constant D (obviously not depending on e). Then we add these inequalities to those of point (iv) of Lemma 
12.41 that we already proved. 

□ 
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2.5 Two standard results 

The preliminaries arc concluded by the statement of two results taken from the existing literature, and that will be used 
in the sequel. 

Lemma 2.7 (Discrete Gronwall). Consider two nonnegative sequences {vn, i»„)„£n, satisfying for some positive constants 

c >o,e e]o, i[, 



Then, for all n [ 



Vn e N*, v„ < v„-i +0v„ +w„. 



n— 1 n 



with Xg := 6/(1 — 0). If w„ = C is constant, we have 



v-a < e"^» 



C 



Proof. Notice first that < e ^ by the convexity inequality of tlic exponential. We hence have 

1 — 6 



V„ < e "Vn-i + ; 



By straightforward induction, we get 



- (1 - e)" (1 - 6i)"-fe+i - (1 - e)" 

ic the following 



n ^ 
-0 + ^(^3^""= 



In case of constant Wn , we use the following fact 
1 



(1 - 0)" 



-C 



1 



(1 - 9)" 



C 



vo ■ 



vo ■ 



C 



□ 



Tile following Theorem can be found in \14f (in the more general case of an infinite dimensional Banach space) where 
it is presented as the "Leray-Schauder Theorem" (p. 286, Theorem 11.6): 

Theorem 2.8. Let (E, || ■ ||) be a Euclidian vector space and 

T : [0, 1] X E — > E a continuous function satisfying T(0, ■) = 0. Suppose furthermore the existence of R > such that for 
any s £ [0, 1], the following a priori estimate holds for the fixed points ofT(s,-): 

T{s,x) = a: => \\x\\ < R. 

Then T{1, ■) : E — ^ E has at least one fixed point in B{0, R). 



3 Approximate system of finite dimension 
3.1 Notations 

We start with a definition related to the discretization w.r.t. time, as in \7t : 

N 

Definition 3.1. We decompose the time interval, ]0,T] = [J ](fe — l)r, fer], where N £N* and r := T/N, and introduce 

k=l 

the finite difference operator : drU^ := . 



We also introduce new definitions related to the reaction terms: 

Definition 3.2. Let us denote by Rf:^ and R^j^ the positive/negative parts of the source terms, and use the same notation 
for . To write things more precisely, we have 

\0 r2j \U2J \0 U2j \S22{U2) + S2l(ui) J 

and obviously 

R{u) = R+(u) - R~(u). 
We introduce the following approximation for this reaction term 

R'^iu) = R+{u) - R-'^u), 

with 



R-'^u) 



Ul \ /7e(sii(mi) + S12(«2)) 



U2J \"/s{s22(u2) + S2liui)) 



where ■y^ is the truncation function that we used in Definition \2. 5\ 
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Let us consider a sequence (Vn)neN of subspaces ofV^{il), such that for all n, Vn is n-dimensional, Y„ C V„+i, and 
UngN^" is dense in L'^(fl). We assume that Vi = (lf2)E, so that the constant function In lies in all subspaces Vn- 

For the sake of clarity, let us denote in boldface the two component vectors: 

\X2j V"'2/ V-^2l("2'"l) 

As previously, all linear operators (such as V,Ot • ■ • j have to be understood line by line in the previous expressions. 
Any dot product (■, ■)-e on some space E of functions defined on Q has to be understood as {u, x)e = (""ii Xi)e + (f 2, X2)e- 
The nonncgative symmetric matrix A'^ (mi , M2 ) will be denoted A^ (it) and we also will often denote by Q'^ (it) the associated 
quadratic form. Finally we will use the obvious notation 



^^(-)=y(n,)j 



in such a way that according to the previous notations, w = 95e(it) and ip^ ^(i") = it- 

We are going to define recursively solutions to a discrete-time problem, all lying in Y„ (n is fixed for the moment). 
Recall that the initial condition it" belongs to L2(Q)2. Given it*~^ 6 L^{Q)-^, we study the following problem (for fixed 
£ > 0, cr G [0, 1], k>l): 

Problem P%(k,u'^'^) : 

Find 11)'° G such that [denoting it'' = ip'^^(w'')] wc have for all x S V^, 

Since E Vn C L°°(f2), it'° takes its values in some compact set of ]0,+oo[2, hence the coefficients of A'^(ii'°) aJJ 
belong to L°°(Q), and tliis ensures that all the previous brackets are well-de6ncd. Furthermore, clearly belongs to 
L°°(r2)2 c L2(f!)2j so tiiat tiie previous inductive definition is consistant. 

The collection of problems (P^ C^^i '"'°~^))fe=i...iv hence a discrete-in-time version of the Vn-weak form of ^ in which 
an e-perturbation and a parameter a have been added. 

We end up this subsection by introducing (bearing in mind the notations of Definition 12. 51 ) tfie; 
Definition 3.3. We introduce the entropy (for some vector function u : ^ ]0,+oopj; 



2 f 

~1 



3.2 A priori estimate 

Proposition 3.4. Under the assumptions of Theorem \1.S\ and keeping in mind the definitions of Subsection \3. 1\ (and 
Definition \2.5\} . there exists a constant Kt depending only on T and the data of the equation (aij , ri, Sij) such that for 
any r > 0, e > (small enough), and a 6 [0, 1] and any sequence (of length less than N) {w^)i<j<k, with solving 
PJ(j, it^^^), the following entropy inequality holds: 



a^,[u'') + raY;^ f (it^KVio^) dx + ^ Hit,^' 



<Kt{<^,{uO) + 1). (14) 



Proof. Plug lu^ in P^(j, it^ ^) to obtain 

of w^-dru^dx + af Q^u^)(Vwi)dx + e\\w^\'li,„= a f kjJ ■ R=(itJ) dx. 
Jn Jn ^ Jn 

Using (i) of Lcmma l2.6l one gets, since Wi = ipf'{ui), 

f ■ d-rU^ dx=- j ■ [u^ - It^-l] dx > -{S{U^) - <?(ltJ-l)). 

Jn T Jn T 



In 

For the reaction term, we have 



/ 10-' ■ i?''(if'')dx = / w{Rl2{u{,ui)dx-\- / wIR2i(uI,u{) 
Jn Jn Jn 



dx. 



For the sake of clarity, wc avoid writing the superscript j for a few lines. Let us focus on the first term of the right-hand 
side (the second one will be similar): 



/ WiRl2{ui,U2)dx = Vl'("l)"l(»'l - 7e[«ll(«l) + ■512(m2)]) dx 

Jn Jn 

= ri ipl' {ui)uidx - / ■!/'i'(Mi)iti7e[sii(Mi) + 512(112)] 
Jn Jn 
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From (mi) of Lemma l2.6l we get the existence of a constant D such that 

xtPl'ix) < D{l + e)(l + Vi(x)), 

which implies 

n / il>l'{ui)uidx<riD{l + e){tM{fl)+S'eiu)). 
Jn 

For the two other terms, we use the fact that ipf ' ^ on [1, +00 [, so that 

- / Vl'("l)'"l7e[«ll("l) + ■512(H2)] dx < - / t/^f'(ni)«l7e[sil(«l) + S12(W2)], 

.'f! Jui<l 

-[ i/)i'(ni)ui7e[sii(ui) + si2(tt2)] dx < (-B + i) / sii{ui)dx + {-B + -) [ Si2(u2)da:. 
in e in e in 

Finally, we use Sij(u) < K{1 + u) to get 

-J Vi'("i)«i7e['5li(«i) + si2(«2)] dx < C ^Aj(r2) + y [«l+«2]^- 
By lemma [2]6] (iv), wo have 

/ [«i + U2] < D{1 + e)ipin) + 
in 

Finally we have, recovering the superscripts that wc omitted before: 

^{£,{u^)-S;{u^-'')) + a J Q^u^)i\7iv^)dx + e\\w^\\l^^^^ < aK{l + ,g;{u^)), (15) 

for some constant K independent of j,n,e,T,cr. Hence if one takes r small enough (such that 1 — tK > 1/2), one has in 
particular : 

which implies by Lemma l2.7l with 9 = C := tK, 

S,{u^) < e^^" [S;{u°) + 1] < e'^-'^(S,{u°) + 1) < e^^^ (<?, (u" ) + 1). 
Plugging this last inequality in II15II . wc get for some constant Ct > 0, 

(T{S^{ui)- S,{ui-^))+Ta j Q=(ttJ)(Vii;J)dx + re||wJ||2i^j^ < Ctt{1 + ^.{u'^)), 
which, after summation over j £ {1, . . . , fc} and using kr < T, ends the Proof of Proposition 13.41 □ 

3.3 Existence and estimates 

The a priori estimate proven in the previous subsection leads to the following proposition about existence: 

Proposition 3.5. Consider the assumptions of Theorem \1.2\ and recall the Definitions of Subsection [37l\ ( and Definition 

For fixed r(= T/N) (small enough), n S N*, to" S and e > 0, there exists a sequence (ti'*)i<fc<]v such that 
solves 'P\{k,u''~^) for all k 6 {1,...,JV}, and which furthermore satisfies the estimate (with Kt depending only on T 
and the data of the equation (aij , ri, Sij)): 

^,{u'') + tJ2[ Q^u''){Vw'')dx + eTj2\\'>"'fHiiS2)^^T{'S'e{u°) + l), (16) 

ll-'llii(n)< [ll-°ll£i(n,+0(v^)]e'-^^ (17) 

fc 

r^||Ji---(u^)||^i(„j < [||^0||^i^^^+O(Vi)] (l+Tre--^^), (18) 

where the sequence of positive scalars {r'^)T goes to r = max(ri,r2) when r — > 0. 
Proof. Step 1: Estimates 

Let us notice that the two last estimates II17I I- II18II are essentially a discrete version of the classical formal estimate 
obtained by integrating our system on £7. In order to recover this rigorously, we notice that since In £ Vn, it is an 
admissible test function for the problem P\{j, u^~^) for all i £ [1, fc], fe £ [1, Af|, and we hence get 



It follows easily from the definition of that the right-hand side of the previous inequality is not larger than Tr\\u^ || j^i ^^^j . 
Thus, 
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We notice then that, 
This leads to 

which implies by Lcmma l2.7l 



ll"'llLi(n) ^ 



k=l 



exp < JT 



1 — rr 



Using the entropy estimate, one has 



Therefore, we obtain 

II'""' III 1 



er^||i<,^||2,j^j </<y(^,(uO) + l) and jT\n\ < \n\T. 



L\SI) ^ [l|-""llLi(n) + ^(^|f^l+^TK.{t«°) + l))]e'- 



with r(T) having the mentioned property, that is precisely JTTJ. 

Quite similarly, wo sum up in II19II and obtain after using the bound on HiwHi: 

fc fe 
Ih'llLMn) +rJ2 = ll^^(^')llLi{n) + ^ (^1^1 + ■^t('?.(«°) + 1)) , 

and we conclude the estimate using the previous bound. 
Step 2: Existence 

In the rest of the Proof, we will work in the finite dimensional Hilbert space E = with the dot product (■, ■)Hi(n) 
and the associated norm. We will proceed by induction and only prove the first iteration it" 6 L^(r2)^ => 3w^ solving 
P^(l,ii''), all the other induction steps will be similar. 

First notice that the problem PJ(l,tt") can be seen in the following way 



Find GE such that Vx 6 E, o-L„o.^i (x) = (x> ■"'"')hi 



where, for in, x £ E, L„o ^ is defined by 



As noticed before, for n; 6 E C L°°{C)^, ip^^{w) takes its values in some compact set of ]0, +oop, so that the coefficients 
oi A^ip-'^{w)) lie in h°°{n), whence L„o^^i £ V*. 

Let us now define a map 

r : [0, 1] X E — s> E 
{a, v) I — > w, 

such that: 

Vx e E o-L„o^^(x) = -e(x.iJ'>Hi(n) ■ 
Such a map is well-defined because of the usual representation theorem for finite dimensional spaces. 
Lemma 3.6. The map T is continuous. 

Proof. Let G E. We want to prove that T is continuous at . Since all norms are equivalent on E, one can use the 
metric | ■ | + || • IIl°°(q) on [0, 1] X E and || • ||Hi(n) E. In all what follows, we assume that v'^ belongs to the following 
open set of E 

\\v'^ - t^^lloo + ||f^ - K^llijl < 1- 

Not that this means that belongs to a compact set of E (just consider the inequalities as large) that we denote K. In 
fact for llxlli/i ^ li oJic has 

\[Ko,^i -I^„o,^2](x)| < ^\\v7^i'"l)-'fi~^{v2)\\i^2^^^ 

+ \\R^v-\v,)) ~ K^(vr'(^2))|lL2(n), 
and one can write for the second term of the right-hand side: 

\\A^^-\v^))Vv^ - A-{^-\v2))Vv2\\^2^^ < \\{A^^-\vi)) - A-(v7i(«2))}V«i|iL2s„ 

+ \\A^<fi-\v2)){Vvi - \'V2}\\-l^2^y 
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Since fj^, A'^ o ipj^ and o ip~^ arc all continuous, they arc uniformly continuous (and bounded) on every compact 
of R^, and particularly on K. Let uik be a continuity modulus for all the previous functions on this compact. Hence, if 
■ui, V2 take their values in K and \\v-i — '^2!l£°^(f2) S: have 

+ ll-4''(</37^(r'2))|lL=°(n)li"i ~ ■"2|lHi{n) 
where we used the continuous injection L°°(r2) > L^(r2). Since all norms are equivalent on E, 

](X)| < CK,n{^K{ri)+v), 

hence using the very definition of T : 

\\T(cr,vi) - T{a,V2)\\-iiUm = - sup [L„o -L„o^J(x) < -CK.n('^x(»?) + »?), 

^ llxli„i<il - I e 

which gives the continuity w.r.t. the second variable. The continuity w.r.t. both variables is then straightforward. □ 

It is now time to use Proposition 13.41 to get an a priori estimate on any fixed point of T{a, ■), for any cr 6 [0, 1]. In 
fact, the case A: = 1 of this Proposition exactly tells us that all this fixed points are in the ball of center and radius 
A'T(f?(it'') + l)/eT. Since clearly T(0, ■) = 0, we now can apply Theorem l2.8l to see that T(l, ■) has a fixed point, which is 
exactly the existence of lo^ and hence the first step of our induction machinery. Inequality I I16I I is then a direct consequence 
of lUll (with o- = 1). □ 

The previous Proposition shows that (for fixed t(= T/N) small enough, n £ N*, -it" £ L^(Sl)^, and e > 0) there exists 
a sequence (n''°)i<fc<]v such that, for all k £ {1, . . . , A^}, denoting := ip^^{w^), we have for all x 6 V^; 

that we also can write (i j = 1, 2j 

drFuu'i: - PnA[af,(«*) + u'[a,j{u^) + eu'^u'^] + ewf - eAiu* = P„ [ijf (n* , )] , 
where P„ is the -orthonormal projection on Vn- 

4 Asymptotic with respect to n 

This Section is devoted to the passage to the limit in the space discretization (n — > +oo ), which can be summarized by the 
following Proposition: 

Proposition 4.1. Consider the assumptions of Theorem \1.2\ and recall the notations of Definitions l275[\3. 31 13.21 and \3. 1\ 

For fixed r(= T/N) (small enough), and e > 0, there exists a sequence (it'^)i<fc<]v of for some p > 1 (and a 

corresponding sequence 

{■w'^)i<k<N) such that for all x G (and for all fc G {1, . . ■,N}), 

^(X«,«^>£2(fjj - ^(x^,u'l'^)l2^^^^- (Ax^, [af, («* ) + u*^ (n* ) + eM^«*] > ^2 (21) 

and which furthermore satisfies the estimates (with Kt depending only on T and the data of the equation (aij, ri, Sij)): 

k 

S,{'a'') + erY,Wf^,^^-^^<KTmu°) + l); (22) 



fc 

Y^j \^M->^'fdx<KT{<Se(u°) + l), 



(23) 



where fia ix) = x 



ll-'=llii(„)<ll-°llii(„)e^-; (24) 

fc k 

i=i j=i 

Proof. We first recall the bounds that hold on it* (and w'^), whose existence is given by Proposition 13.51 
Using lfT2t . for all fc 6 {1, . . . , N} and n S N*, we have 

Q£/ fc\/r-7 fc\ \ ^ £ ff ^,rL\\r-! k.n\2 , ^ £ '/ k,n\\r-! k,ni2 
KJ(Vii^„) > -^a|i («i' )|Vui' I +-j^af2(«2 )|V"2 I ■ 



Since afj(a;) = aij{x) + ex, assumption H2 leads to 

Q^u'^){Vw'^) > — ^ I V/3„(«f ■")|' + I V^.(n^'")|' + 4£ 

(1 — a)^ I I (1 — a)'' I I 



V/ fc ,ri 
V 



2 



■Ae 



V/ fc.n 
V"2 



(26) 



11 



(where we recall that I3a{x) := x 2 ). We also have, because of I I17I I. the boundedncss of (ttjij)„gf!j« in L^(f2). Since 
the asymptotics that wc arc studying is only w.r.t. n (that is, e and r are fixed), we see finally that j y u 



and 



nSN* 



^2 



are bounded in H-'^(Sl) ^ ti^ (C), with 2* = 2d/{d — 2) > 2, and hence (ttJ^)„gN* is eventually bounded in 



some \J'{Q,) space with 1 < p < 00. 

On the other hand, estimate II13II gives us 



V"2l(«J'")'^f2(«2'") 



vhich together with III6I I leads to 



1 



VV''J2l("l'")»12("2'") 
Using Poincare-Wirtingcr inequality, we get, for some constant Cn depending only on 

< Cn 



/ =- / k.n\ fr I k,n\ 
V'*2l(«l' )"12("2 ) 



V'*2l(«l )«!2{«2 ) 



< C'n [ll«2i("i'")llLi(n) + ll'Ji2(«2'" 
and, again because point {iv) of Lemma l2.6l for some constant D {e < 1), 



^ KT{<ge(uO) + I) 



< Cnyila2l(«t'")llLl(n)IW!2(«2'")llLl(f2) + ^^\/^T{'^e{uO) + 1) 



2v^ 



^^iCT(<ft(«0) + l), 



so wc finally have 



L^{n) ^ 2\/t V 



and hence thanks to III6I I. the sequence ( \/ a|-^ (it-^ '")aj2 (^'2'") 



20= 

is bounded in H^(f2) 



L2*(n), with 2* = 2d/(d- 



2) > 2. The previous continuous injection implies then that (alj («j '")af 2 (^^2 "))"6N* is bounded in some L*'(f2) space 



with 1 < p < 00 and since iji('i*''") = o,ji{u^'") + eu^'", we eventually get that (wj'"ii2'")neri 



,"1 ai2(«2 ) 



1 / k.n I fc,n\\ 

and I 'J2i('i'i ) I 



are bounded in the same L*'(f2). 



Summing the bounds already obtained [we recall that they hold for a given e and r], we see that (for all k G {1, . . . , A'^}): 



k.n 



and 



k.n 
'^2' 



/«2l("l'")«12("2'") 



arc bounded in H^(f!), 



is bounded in H^(f2), 



/fc,n\ I k,n\ I k,n k,n\ 

(.f 1 jnSN* I (.^2 jnSN* i 1*2 ^lEI* 



ai2(«2 ) 



■1 / k.n I k.n\ 

and (121 ) 



L*'(f2) space, with 1 < p < oo, 

and obviously (io*)„gf!}» is bounded in H-'^(f2), because of estimate III6I I. 



n6N* 



are bounded in some 



Since we are only dealing with a finite number of values for k 6 {1, ■ ■ ■ , A'^} (at this point, the functions are not 
time-depending), we shall only detail the study of (it^)„gK*i the other values of k being similar (one only has to extract 
a finite number of subsequences). For every test function x = (xi)X2) 6 V^j, the weak formulation II2UI I may be written 
(i^i6{l,2}) 



1/ 1/ 0,7i\ /A fe/liTi-Ni 1.^^ /l<7i\i ^■^n \,rL\\ 



(n) 



First we extract (but do not change the indexes) a subsequence of (io^)„gN* converging in h'^(fl) and almost everywhere to 
some element S H^(f2)^. The almost everywhere convergence is transmitted to (u^)„gf^» , since ii^ = ip^^{w^), the ipf 
functions being homeomorphisms. Because of the L*'(f2) (1 < p < 00) bound for (ii^)„gKi*, we can extract a subsequence 
(of the previous subsequence) converging weakly in hP{Q) to some function 6 L^(0)^. It follows by a classical argument 
that is almost everywhere equal to := (p^^{w^), that is the almost everywhere limit of (itJj)„gN*. The same 

argument holds for (ji]^'"u2'")„gN* and (u^'" aij(u^.'")j , i 7^ j = 1,2, since we have the same LP(n) bound. We may 

also assum (Qthat (MO'")„gN- converges weakly in YJ'iVC) to some function tt" £ L*'(S1). Eventually, the cutoff perturbation 

introduced on and the (superlinoar) reaction terms ensures the weak convergence, in L^(f2), of (ifi("i'")) and 



^In fact we have u"'" = u*^* \ This line is just to justify the handling of the same term in the other time steps. 
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{R"^ {u^'"))„£m*) respectively to cif^{uj) and R'^{u^). 

We now fix a test function x in U„gN* V^, say x S V^, m E N*. Then, for n greater than m, since the sequence of spaces 
Vn is increasing, we have the previous weak formulation that passes to the limit thanks to all the extractions that we 
made, and get eq. I I21II . 

It remains to prove the bounds announced in the Proposition. 

Since we have almost everywhere convergence of (itJi)n6N a^nd weak H-'^(C) convergence of {i«J^)neNi Fatou's Lemma 
and the classical estimate for the weak limits give directly, for all k S |1 , , estimate II22II . 
We use then II26I I that we rewrite here: 

Q^u'^)iVw'^) > —±— + |v/3„(«^-'")|' + 4e 

(1 — a)-' I I (1 — a)^ I I 

where we recall that I3a{x) = x 2 , We notice that for all j S the sequence (■i'n))neN was, precisely in view 

of the previous inequality, bounded in H^(f!), so that we can (adding another extraction) assume that we also had the 
convergence (/3ct (tt^))neN ^ I3a(u'^) in H^(f2) (using the uniqueness of the weak limit) and we hence have estimate II23II 
using Fatou's Lemma. 

As for III7I I and III8II . we notice that we had a L''(f!) (p > 1) bound for both {u^-")„^m* and {R~''^ {u^'"))„eM* , 
associated with almost everywhere convergence. This is sufficient (using EgorofT's theorem) to get the (strong) convergence 
of these two sequences in L^(Q). Estimates 1178 and Ill8t give hence Il24t . and Il25t after taking the limit in n. □ 



'4e 



5 Duality Estimate 
5.1 Notations 

Definition 5.1. For a given family h := (h^)i<k<N of functions defined on fl, we denote by hJ the step (in time) 
function defined on ]0,T] X Q by 

N 

}f{t,x) := X]'^*(^)l](fe-i)r,fcr]W- 

k = l 

We then have by definition, for allp,q £ [l,oo[, 



1/9 



and in particular 

i/p 



5.2 Duality estimate: abstract result 

This subsection is devoted to the cstabUshmcnt of a discrctizcd version of the duality estimates devised for singular paraboUc 
equations in \23f . \21\ . We start with the following Lemma: 

Lemma 5.2. Consider a smooth function b S 'i!f°°(f2) such that b > 7 > for some constant 7. For a given ^ £ L^(r2), 
the variational problem associated to the equation 

-•!> + bA<S> = (27) 

is well-posed in H^(Sl) (with L'^(Sl) test functions) . If 'I' is assumed to be nonpositive, then the solution of the previous 
equation is nonnegative. 

Proof. Since b is bounded from below (by 7) and above, we can first (uniquely) solve the (variationnal formulation of the) 
problem 

-fe"^<I> + A$ = 6"^*, (28) 

in the Hilbert space H^(r2), using Lax-Milgram Theorem. 

Elliptic regularity ensures that the constructed solution lies in fact in H^(Q), and even in H^(f2) because of the 
variationnal formulation, so that l(28j is satisfied almost everywhere and we recover a solution of II27II after a multiplication 
by b. Uniqueness is straightforward. All L^(r2) test functions are obtained by density. The last part of the Lemma is 
obtained by maximum principle. □ 

We then turn to the: 

Lemma 5.3. Consider a real number r > such that 1 — 2rr > and two families b := (b'°)i<fc<]v, F:= (F'=)i<fc<]v of 
functions. Assume that for all k e [1, 7V|, b^ >l and F'' < (pointwise). 
Then there exists a family $ := ('J'''")i<fe<jv S of nonnegative functions such that, (defining ^^'^^ := 0), 

Vkell,Nl ^ ~ +b''A'S>'' = Vb^F''-r'i>'', (29) 
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where eq. 129\) has to he understood weakly, against test functions. 

This family satisfies furthermore 

VjG[l,7V], l|V$^||^.(,,)<e'-M^||F-||^.(^^,, (30) 
||vrA*-||^.^^^,<e'-M-||F-||^.^^^,, (31) 
Vie I1,JV], ||*^||Hi(n) <e'-M^C^o[e'-W^+||v^||i2(Q^)]||£-||i.(Q^), (32) 
where the sequence of positive scalars ('''(r))r goes to r when r — > 0. 

Proof. A straightforward induction using Lcmma lS. 21 gives us step by step the existence of the family {^'')i<k<N a^nd the 
nonnegativity of its elements. Now in the fc-th equation, rA^*"" £ L^(f2) is an admissible test function and we have hence 
after integration by parts: 

/ |V$'=pdx- / V<I>''+^ ■ V-I-'^dx + r / b'=|A$'=pdx = r / vV'F*^ A^'^dz + rr / |V#'=pdx. 
Jn Jn Jn Jn Jn 

Using 



and 



we get, for all fee [1, A''] 

/ IV^^pdx- / |V*'=+l|2dz + r / b*|A*'=|2dx < r / \F^\'^dx + 2rT [ IVJ^I^dx. 
Jn Jn Jn Jn Jn 



In Jn Jn 

Let us introduce the auxiliary sequences (1 < < A^): 



(l-2rT)'=/2' 



• (1 _ 2rT)'=/2 ' 
we have then, dividing the previous inequality by (1 — 2tt)^'^^ , 

i / |V*'=pdi:- / |V*'=+M2dx + / 6''|A*'=|2dx < / \G^\'^dx 

l~2rT Jn Jn l-2rT Jn - l~2rT Jn 

1 - 2rr in 

hence 

/ |V*'=pda:- / |V**+^|2dx + / b'^IA't^l^dx < / IC'^l^dx. 

Jn Jn 1 — 2rr Jn 1 - 2rr Jn 

Now for J 6 [1, N\, summing up over j < k < N and using $^+^ = 0, we get 

/ |v*^fdx + -^^ / 6'^|AVI/'^|2dx< / IG^l^dx, 

Jn l-2rr^Jn l-2rr^yn 

and hence, in terms of the 

which gives immediately I I30II and JSTJ in the particular case j = 1. 

Now for the last bound, notice that after integrating the fc-th equation, we have 

(1-rr) / <I>''dx = / $'=+Mz + r / A** - V^F''} dx. 
Jn Jn Jn 

Defining this time 



(1 - rr)'-- ' 
pk 

(1 — rr)*" 
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we get, dividing the previous equality by (1 — rr)'''^^ , 

I ^''&x= f **^+MxH — / {6*^A^'''' - V^G'=}da;, 

Jn Jn 1 - rr Jq 

which, after summing up all over j < k < N, leads to 

N 

Jn ~ Jn 1 - rr in 

1 ^ /■ ^ 

k — j 

and we have hence, using (131t . 

/^*^d.|<[||v^||L.,^^,+e^M-]||£-||L=,^^,e^M-, (33) 
so that I I32II follows from 1301 1 and Poincarc-Wirtinger inequality. □ 

If we want to use functions of the previous family $ as test functions in the weak formulation II21II , we have to show 
that they actually belong to W^'^ (Q). This can be done using the following lemma 

Lemma 5.4. Let r be a strictly positive real number. If w £ Hj(f2) satisfies w > and —Aw < f + rw almost everywhere, 
for some f in L''(H), q > max((i/2, 2), then 

\M\L°^(n) ^ C'n,q,r (ll/l!£9(n) + Ihli ii(n)) • 

Proof. Cfj.q and Cq denotes constants that may vary from line to line. Lemma 15.21 gives us the existence of g £ H^(f2), 
unique solution of 



Ag = {r + l)w + / e L™'"(2V9)(f^)^ 
ng W^-P{n) ^ hP' (n). If 2* < q, th 

\\g\\^2"i^) < Cf2||g||^2,2«(n) < Cn\\(r + l)w + f\\^2' 



where p* is defined by the Sobolev embedding W^'''(f2) ^ (f2). If 2* < q, then the elliptic regularity of (Id — A) and 
a Sobolev embedding (cf. |14p ensures 



But we also have 

w — Aw < (r + l)w + /, 

that we can write 

(w-g)- A{w -g)<0, 
so that by the weak maximum principle in Hj(f2), we have 

< to < 5 e L2**(n). 

The previous argument can be iterated in order to finally get w £ L'(f2), and thus using the elliptic regularity mentioned 
before, 

Il9llw2.9(n) < CnWir + l)w + /||L<!(n)- 

Since q > cZ/2, we get the h°° estimate 

llffllL~{n) < C'n,9ll(^ + l)'" + /llL'(n)' 
that is again transmitted to w, thanks to the weak maximum principle : 

\ML^(n) ^ C'n,g||/|lL9(n) + Cn,q{r + l)\\w\\^g^^y 

Since 

II iiL''{n) - II iiL^(f2)" "L=°{n)' 

we get by Young's inequality 

Cn,9(f + l)ll«'llL«(n) < ^'lAMli.^n) + ^ll"'llL'=(r!)' 
and we are eventually able to conclude that 

ll«'llL°'(n) < (ll/llL«(n) + ll'"llLi(n)) • 

□ 

We now can prove the following: 
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Lemma 5.5. Under the assumptions of Lemma \5.3l the sequence of functions ^ satisfies in fact, for all k 6 [1, A^l 

||<I.'=||i»(„) + ||A$'=||i^(„)<C.,a[l + ||v^||i2(Q^)]||£-||i^(Q^), 

the constant C-r,n depending only on r and Q, but (severely) blowing up as r 0. In particular we get that £ W'^'^(n) 
for all q a [1, oo[. 

Proof. Since (I-^ > 0, = 0, < and b'^ > 1, the Af-th equation of l|29ll gives us 

and since S L°°(f2) > L*(f!) for all s > 1, the previous Lemma applies: 



ll*^llL-(0)<Ca4||Zl|L~(Q,) + ll*^llLi 



and we use then II33II to get 

ll-f'^llL-to) < C.,n.T[l + ||v^|Il^«5^)] 
For all k S [1 , A''] , using the fc-th equation of II29I I , we have in the same way 

and hence a descending induction gives eventually (C,- n varies from line to line) 

ll*'=llL-(n) < Cr,n [l + ||vniL-(Q^,] ||£1Il-(q,)- 
To handle the laplacian terms, we write for all fc £ [1, A^l, using again II29II . 

(J)fc+i _ <j)fc i?* <j>fc 
-A*'' = ■ — + r — , 

which directly gives, using the previous estimate, 

W^^'Wl.-^n) < Cr,n [l + II v^|lL^(g,)] l|i^1lL~(Q,)- 
In particular, the norms |j<I>* — A<I>'^||]^gjQj are all finite for q £]l,oo[. Thereby, using again the aforementioned (see 1141 ) 
result on elliptic regularity, we get the finiteness of all the norms Ii'l'''|iw2'!(f2) all g < oo. □ 

5.3 Duality estimate: application to the system 

Proposition 5.6. Under the assumptions of Theorem and using the Definitions \2.5\ \3.1\ and \5.1\ any weak solution 
of eq. 121\} satisfies the following estimate: 



< D, (34) 

L\Qt) 

where D is a constant depending only on T, Q and u° . The same estimate holds when the subscripts 1 and 2 are exchanged. 

Proof. We shall denote by and n 2 the step (in time) functions (see Definition 15.111 associated to the two families 
('^i)i<fc<]V a-nd ('*2)i<fc<iV (excluding hence and M2). Since the original system is symmetric in ui and U2, let us focus 
on Ml for now. As in Lemma 15.31 we consider two families b := (fc'')i<fc<]v and F := {F'')i<k<N of "^""(H) functions 
such that 6*° > 1 and F*^ < 0, and the associated sequence <£> := {^'')i<k<N 1 with r := ri. As shown in Lemma 15.51 the 
functions <E>''' belong to W^''' (f2) (for all p > 1) and are hence admissible in the weak formulation I I21I1 . We may therefore 
write, for all fc 6 |1, A'], taking for test function in the weak formulation on mJ: 



L^(n) 



Thanks to the regularity of one can take mJ" as a test function in the fc-th equation II29I I. so that 

-?)L^a) - ;<*''«^)l^(0) + <^*''f''-^)L^O) = (^^'.«?)L^(n) -^)l=(«)- 
We hence have, for all fc £ [1, A'], adding the two previous equations. 



recalling that 0^2(2;) := ai2{x) + ex, a\^{x) := xd\-^^{x), d\-^ := 7e(dii), and the decomposition = -R12 ~ ^12^ ■ Since 
is nonnegative, we get, if we denote c^ := 6* — cZjj('uJ) — a|2(«2)i a-^er summing over k £ [1, N\ (and recalling that 
$iv+i =0), 

. JV N N 

fc=i fc=i fc=i 

The following approximation Lemma will help us to handle the sequence c := (c''")i<fc<]v: 
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Lemma 5.7. There exists a sequence of families (6m)meN ~ (^mi ^mi • • ■ > ''m)™sN of'^°°{Q) functions such that 

V(fc,m)G [l,iVl xN, bt>l, (36) 

0, (37) 



VfcG[l,iVl, / K,-dUu'i)-aUui^)\(l + ut)dx 



m^-\-oo 



(l + u'l)dx — y 0, 



in 

/ \dli{ui) + al2{u2)\dx = 



7?l— > + 00 



I L (Qt) + = 



(38) 
(39) 



Proof. This Lemma is a rather direct consequence of the regularity of the boreUan measure we know that 

cZ| (Uj) + 0^2(1*2 ) S L^(0, (l + nj)da:) and is lower bounded by a strictly positive constant, we may hence approximate it in 
this space by smooth functions, which are still lower bounded by a strictly positive constant. This gives II36I I- II37II . 

Using the mentioned lower bound, we have 

i/f4-^df,(«5;) + af2(n|) < " /^LKH^f^^ + V'^fi K) + «!2 («2 )] 

= \bi^-d\M)-'^l2{y-2)l 

so that we also get II38I I. Noticing that 

II I — ii2 ^ r 

we see that II39II is a simple consequence of the convergence II37II . □ 



We now fix a family F := (F*-') 

l<k<N of nonpositive functions, and for each family bm of the sequence 

{bm)meN defined in Lemma 15.71 we define the corresponding family $m := (^^^j, <I>^, . . . , <E>^), using Lemma 15.31 The 
previous estimate II35I I can now be written 



JV N N 

k = l fc = l 1 

where := 6jjj — d^-^(tij) — a^2('^2)' Since the right-hand side is nonpositive, we may write 







(f2) 



Using II32I I of Lemma [53] and II39I I. we get (for a given r, e). 



L^{n) 



+ \\\ k 



l|£llL^,g.)^'-^<^'" 



— > Cn\\u lii 2,„, 

m->-|-cxD {'•') 

Using Lemma 15. 51 we get on the other hand 



='-i(t)T . 



/dl^{ul)+al2{ul) 



||T|| < T7V^^sup^{||A<I.^|lL^(j^)||c^nJ||j^i^^j] 



1+ \/&' 



Iz-I 



l°°{Qt) 1 II'''"''! llLi(n) 



using II37II and II39II for the convergence. In order to treat | 4 | , notice that wc know, from inequality I I22II : 
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hence, using again I I32II of Lemma l5.3l we get {tN = T) 



e'-i(-)T + 






















e'-i(-)T + 











(n) 



|£llL~(Q,)-=v^{Eli"'?liHi(0)} 



1/2 



l^"llL°°(Qr) 



£Tji^T('S'eK) + l), 



and we may use again I I39II to get that the last upper bound is converging, as m goes to +00, to 



V4i(«T)+'^12(m^) 



Finally, because of II38I I. wo can see that | 1 | converges, as m goes to +00, to 



N . . 

/ u'i^dl^{u1)+al^{u!l)F''dx = / ^dfi + af ^ (^^ )£-dx di 



All the previous estimates give hence, denoting ^Tv'^ii (l^ii) +'^i2(— 2)' 









JQt 







X [V^^KriS^iuO) + 1) + ||w°||l2(„J ||£^||l= 



(Qt)' 



Since ft'^ is a step (in time) nonnegative function and the previous holds true for all non-positive smooth (in x) step (in 
time) functions, we have then by duality 



< e 



2ri(r)T 



Cn 



1 + 



X [x/irv/i^T(<f.(n«) + l) + J. 



Hence, in norms, we get 



Cn 



l + IHiibT) +a!2 (2^1)11^1(0^) 



Since e: < 1 and {ri(T))r — > ri, the previous inequality can be written in a simpler form : 



< D 



(Qt) 



1+ ||rfll(MT)||Ll(0) + II«12(M2)|Il1(Q^) 



where the constant D only depends on T, SI and it". Because of the concavity of (and since it vanishes in 0), we have, 
with some constant Ci,C2 not depending on e, 

af2(x) < Ci +02 1/^1 (x). 



so that, from I l22t . we have 

kl^ +ai2(W:2) T 1 

ML J lllj (Qt) 

that we may write (changing the definition of D) 

||«T^[dii(«I) +ai2(M2)]||j^ij^_^j < D 



1 + ||dn(lti)||£^i 



(Q) 



18 



But wc have then, since d?, < dn which is a non-decreasing function, 



||2iT^[dii(ltI) +al2(2il)]j|j^i < -D 



< D 



1 + |i2i->v^'^ii(^iT)|f_i,„^ + \\'^uri<V2Ddli{u'i)\\^i 
1 



lL^{n) 
lL\n) 



L^(n) 



so that 



mT'[J'^!i(«T)+»!2(m2)1 , <D\i + di^(V2D)fi{n)], 



that is to say, eq. 1341 1. □ 

6 Proof of the main Theorem 

This section is dedicated to the proof of our main Theorem. 

Proof. Thanks to Proposition 14. ll wc can consider a solution u\'^,ul^'^ to system II21I I. For the sake of clarity, we shall 
denote this solution by tt^, 2i 2 • keep in mind that since we are dealing with the limit (t, e) — ^ (0, 0), we shall have to 

use bounds that are uniform w.r.t. these two parameters. Hence in the sequel, if not mentioned, the term "bounded" will 
always have to be understood as "uniformly w.r.t. r and e". We shall only work with u'^, the study of m 2 being exactly 
identical. 

We start with a simple computation. Take some real number r/ > and write 

JV 

Ul{t + V,^)-U\{t,^) = X! [l](fe-l)T,fcr](* + »?) - l](fe-l)r,fcT](*)]"l(a;) 
fc = l 
JV 

= X! [l](fe-l)T-J7,fcT-T,]{*) - l]{fc-l),feT](*)]«l(^)- 

k = l 

Hence, using 

l](fe-l)r-T;,feT-77](*) ~ l](fe-l)T,fcT](*) = l]{fc-l)T-J7,(fc-l)T] (*) — l]fcT-r;,fer](*)i 



we get 



u'[{t + rt,x)-u[{t,x)= ^ l]feT-,,,fer](i)«J^^(a:^) - llfeT-T,,fcT](i)"i(a;) 

fc=0 k = l 

JV-1 

= l]feT-„,fer](t)[«J"^^(a:) - u'[{x)] + li_^,o]W«i(a:) - llT-^,T]W«f (a;)- 



If we denote by c,, the translation operator g{-,x) 1— > g{- + r],x), wc have, introducing the space E,, := h^{[0,T 
r,];W2'°°{n)'), 



1 iiw2-°= (n)' 



Now recall that, for all j e [1, N}, 

weakly We then have (taking j = k + 1), since }1^{Q) ^ R-^{n) ^ W2'°°(n)', and L^{n) ^ W2>°°(n)' 
= l|[rf!lK+^) + af2(4-+i)]«f+il^.+i^ 11 



(n) 



^1 ^illL^n) 



<i|lLi(n) 



+ eCo|hJ+'|lHi(n) +llfiI2"(^^'+')liLl(n)+^•lli"^^'llLl{n)• 



Using In = l^fc+i^^ + l^fc+i^;^, we have 



KiK-+i) + af2(«2*+i)]«J+i| 



iLi(n) 



2 



sup rf!i(x) + ||af2(n2 )llLifn)- 
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Since for e small enough, coincides with dn on [0, 1], using also the last point of Lemma 12.61 to get D(2 + ^2)' 

we have eventually (recalling ((22}) for some constant C not depending on e, 



< 



where we change the definition of C, so that finally 

N-l 



L^(n) 
2 



k = l 
N-l 



-c, 



+ C + £Cn||t«J+'|lHi(n) 



<4r, nW'ifiC"!) +a!2(2i2) L,^ + r,TC + ryerCn ^ ||ti,'=||Hi(n) 
II IIL (Qt) 



1\ — ± 

-ri E -{ll^^2-(-'+^)llLHa) +-ill"i'+'llLi(n)}' 



which, using 1134 l l. II22II . II24II and II25I . is going to as r; — > 0, uniformly w.r.t. t and e. 
We also have 

N 

|v«ii(i,x) = J2 h(k-l)T,kT]{t)\'^u'[\{x). 

k = l 

1-a 

Recall the notation I3a(x) := x 2 , Given v : Q ^ M.-^- such that /3ct(t)) G H (f!) (and considering for example the sequence 
{v + l/m)^gN), we get 

1 — Of Q + l 



weakly. Hence, using II23II . we have for all k 6 Jl, A'']: 

fe 



■J2 f i«ir"-'iv«iP 



dx < Kt{<^e{u°) + 1). 



We know that 



|V«I|(i,a;) = X^l,(fe_i),,fe,,(t)|nJ|^hJr^|V«J| 



Hence, using the dual estimate I lij4t . we see that Vm^ is bounded in L''° (Qt)i where is defined by the equality (a 6]0, 1[) 

0+1 1 _ 1 
4 2 ~ ^' 

that is (Jq = 8/(2a + 6) > 1. At this stage we obtained, denoting := L^{[0,T - ri];W'^'°° {n)') , 

• ||a'^ii'5^ — m^IIe^ goes to with rj uniformly in t and e, 

• {uD-r is bounded in L9° ([0, T]; W^-''" (H)) . 

Wc know that W"'^'''°(f!) ^ L'~(Q) > W^'°°(n)', where the first injection is compact, the second one is continuous, and r 
denotes any real number of the interval [l,q*[. We may hence use Theorem 5 of 1251 to get the compactness of {u\)t in 
Li([0,T];L'-{Q)). 

From the definition of : 

N 

2il(t, a;) := u'l{x)l](k-i)r,kT]{t), 

k = l 

we get 



(T—rU\(t,x) ■.= u\{t-T,x) = EMl'{a;)l](fc_l)^,feT](t-T) 

k = l 
N 

= ^■^li^)hkT,{k + l)T]{t) 

k = l 
N+1 



E mJ ^(a;)l]{fc-l)r,feT](t)- 
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We hence obtain the following expression for the rate of growth 

^ x) = J2 l]{fe-i)r,fer](*)^-"J W + -ho,T]{tH{x) - -l)T,T+r)(*X(^), 

''' fc=i ^ 



recalling the definition of the finite difference operator, 

:= ^ 

Now let us fix a test function 9 e &{[0, T[; V;^{n)). We have 



/ / ^1 " g--rMi = / {d-ru'l,e{t)) At+ - {u\.e{t))At, 

Jo Jn r ^-'{fc-l)r T Jo 

where we denote by (■,■) the duality bracket S^' (Q) / !^{Q) (which is simply the integration on Q here ...). A change of 
variable leads to 

/ / =i ■e{t,x)dxdt= / u[- -{t,x)dxdt- / u[ ■ — {t, x) dx dt 

Jo Jn T Jo Jn "T J~T Jn t 

= / ul —{t,x)dxdt 

Jo Jn T 

m 



/ / ■ — (t, x) dx dt — f f u\ ■ —— (t, x)dx 
Jt-tJq T J-rJn T 



dt 



Then, | 1 | equals to for t small enough because 9 S ^{[0, T[xf2), and | 2 | equals to for all r because of the definition of 
It J, so that finally 

y {d-ru'l,e(t))dt+- r {u\,e{t))dt= r ^ [ u\-^—^{t,x)dxdt. (40) 

^i{fc-l)r T Jo Jo Jn T 

Since for all fc £ [1, N\ we have (in the weak sense) 

dru'i - A([dfi(«f ) + a!2(«2)]«j) +e(m'= - Au-f ) = flfaC'"'"'): 
we also can write, for all t 6 [0, T[, 

{d.ru\,e{t)) - (af 1 {u\) + ai2(«|) + eu\ul, A9(t)) + e{wl,e{t) - Ae{t)) = {Rlzi^''), e{t)), 
so that integrating on ]{k — 1)t, t] and summing over k £ [1, Af], we get 

/ {dTu'l,e{t))dt- [ f [afi(«'^) +«'^ai2(M2) +e«lM2l ■ A9(i,x)dxdt 
^iJ(k-i)T Jo Jn'- J 

+ £V / ,0(t) - A6»(t)) dt = / / R\^(u;'){t,x)e{t,x)dxdt. 
^r^i(fe-i)T Jo Jn 

Using 1401 1. we eventually get 

/ / l[o T— r[(t)lii • '^{t,x)dxdt- / [afi(M'[) +?i'5^ai2(?i2) +£2iTii2] ■ ^^(*'^)da::di 



N .J,,- 

-e^/ («,f,e(t)~ A0(t))dt 

- f {ul,e{t))dt+ f f Rl2iyr)it,x)9(t,x)dxdt. 
T Jo Jo Jn 



We now can study the limit (r, e) — > (0,0). Thanks to the compactness result that we proved above, we get the existence 
of It := («i,'U2) G ([0, T] ; L''(f2)) (with r < <?*) such that, up to a subsequence, we have 

(iil)r,e > Ul, 

(r,e:)^(O,0) 

(liD-r.e > U2, 

(r,e)^(0.0) 

in ([0, T]; L' (f2)) , and also almost everywhere on Qt- Because of the dual estimate 1341 1 {^i^\J 0,12 (1^2)) bounded 
in L^(Qt)- Using assumption H2, we see that ai2 is at most linearly growing, so that using again the dual estimate 13411 
(but inverting the subscripts 1 and 2), ^ai2(M2)) bounded in 1?{Qt)- Writing 



2ila,i2(2i2) = 2ii\/<ii2(2i2)\/ai2(l!:2)' 
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we see eventually that (^^"^0,12 {u^)^ is bounded in L^/''(Qy) and we may thus extract a subsequence converging weakly 
in this space, and whose limit has to be equal to Uiai2{u2) (because of the previous almost everywhere convergence). 



As for the self-diffusion, the dual estimate II34I I ensures that y^dfj"(«Y)') bounded in L^(Qy). Since d^-^ is 

a nondecreasing function, wc infer the boundedness of ^cifi(lt^)^ = ^mT'^ii (^il)) iii ^^{Qt)j so that the weak 

convergence of this sequence can be deduced from its uniform integrability (Dunford- Pettis theorem). Since df-^^ < dn 
which is a nondecreasing continuous function, introducing the pseudo-inverse 311 (t) := inf{x 6 M.^; x dii{x) > t} which is 
also going to +00 with t, we get 

^«T'i!l(«I)lu-jdfj{u-j)>M(*,a;)dxdt < ^ "1^11 („-j)>A/(t, a;)dxdt 

= / / IiIrfii(liI)l<>9ii(A/)(t.a:)dxdt 
Jo Jn 



- gii(M) 

which indeed goes to with M , uniformly in e, r, thanks to I I34I I. 
We hence get 



■ /^e{t,x)dxdt. 



Since d is smooth, wc have 



-dtb 



uniformly on Qt, and 1[q ^_^[ converges to 1 in all L°(Qt) (s < 00), hence 

/ / llOT-rl(t)«I —{t,x)dxdt > -/ / Ul-dte{t,x) 

Jo Jn ' r (t,£)-»(o,o) Jo Jq 



dx dt. 



For the reactions terms, since the nonlinearities are always strictly sublinear or dominated by the self diffusion, one easily 
manages (using the dual estimate) to use Dunford Pettis criterion. Indeed, first, since (7£)e is increasing to the identity, 
the uniform integrability of Ri2^{u^) reduces to check this property for both sii{u'^)u'^ and si2{u'^)u\. Since dn is 

nondecreasing, wc always have lim — ^ \ = 0, so that using again the pseudo-inverse /ii(t) := infix 6 M.^; x sii(x) > 

z— > + oo zdii (z) 

t}, which is also going to -|-oo with t, we get for M > 



/ / ■5li(2iT)MTlsii(u^)u^>M(*,a;)dxdt = / / sii(«I)?iIluT->/ii(A/)(i,a; 
Jo Jn Jo Jn 



< sup 



siiiz) 



>/ii(A/) zdii{z) 



u\Jdii{u{) 



)dxdt 



L"(Qt) 



which goes to with uniformly in e, t thanks to I I34I I. We used here the fact that estimate 1341 1 is still true if one 

replaces d\-^ by d\\ (Fatou's lemma). For si2(l'2)2ii> using again the pseudo-inverse, it is easy to exhibit some positive 
function I going to -|-oo in -|-oo such as 512(^^2)^^! — large enough, < £{M) => > £{M). In the case of the 
previous implication, one may easily use the already noticed uniform integrability of (m'^)t and if ^ ^(^)i then we have 
by Young's inequality for any small <5 > 0, 

Ul si2(«2) <Sul ai2{u2) + ■ *12("2) 



5 012(112) 



so that we may conclude using I I34II and assumption HI 



lim 



^12(2) 



*+°° 2 ^^^22 (2) + 012(2) 



: 0. 



For the other terms, first notice that 



{wl,e{t) - Ae(t)) dt 



fc=i 



JV 

< e\\e - mi^^^yo^TV^.\n)) J2 ^ll'^^llL^n) 

AT 
k = l 



< 



< 



where we used 12211 for the last but one inequality. 



A0| 



1/2 



lie - Ae| 



L°°([o,T]iL^{n)) 



* (r,E)^(0,0) 
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<4yrih\Qjhl\\L\Q^)\\^^h-'iQry 



dx, 



Wo also have 

/ / eu\u'^Ae{t,x)dxdt 
Jo Jsi 

which goes to zero with (r, e) thanks to the dual estimate J34j. 
Finally because of the continuity of 9, we can write 

- r {ul,e{t)) dt = r {ul,e{Tt)) dt — > {u1,e{Tt)) = [ u'l{x)e{o,x) 

r Jo Jo (T,e)^{0,0) 7n 

so that it only remains to show 

- r {u\-u\,e{t))dt= r {dru\,e{t))dt 0. 

T Jo Jo (T,e)^(0,0) 

From equation II21I I for k = 1, i = I, j = 2, we have (xi = S{t)) 

{d^u\ , e{t)) ~ [af 1 («; ) + u\ai2 («2) + e«i«2] . '^^'W) + ^i'^l . ^ (*) " ^^{1)) 

= (/?f2{n;,ni),0{t)), 
that gives, since L^(f2) ^ W^'°°(r2)', for some positive constant Cg depending on 9, 



1 



/ {drul,9{t))dt 
Jo 



< T\\al-f^{u\) + ulai2{ul) +eulul + ewl + -Rf2("i. "2)liLi(n) 

< ||2iIf^!l(liI)l[0,r]xn|lLl(Q^) + il«I'^12(li2)l[0,r]xnllLl(Q^) 

+ 4vri\\l2^Q^^+e\\yr2\\l2^Q^^+eT^/JIia) 

As shown before, all the sequences written in front of the characteristic function l[o,r]xr! ^'''^ uniformly integrable (we 
actually showed this for the self diffusion term, and have some h^(QT) with p > 1 for the two others), and they hence 

vanish with (r, e). The three other terms also vanish, because of the duality bound for the ||«7I|t ■> terms and because 

' ' i^~(.Qt) 

estimate i22t for the I1'''iIIl2jj^j term. 

We get then 

/ ui\dte~[dii{ui) + a22iu2)]Ae]dxdt= f u1(x)9{0, x) dx, 
Jqt ^ ' Jn 

that is the weak formulation on Qt of the equation 

dtui - A|^aii(ui) + «iai2(ti2)] = Ri2{ui,U2), 
initialized with ui{0,x) = u^{x). A similar (symmetric) formulation holds for U2. 

We end up the proof of our Theorem with the passage to the limit in the duality estimates. First notice that the 
constant D in II34II can be written as a polynomial function (with positive coefficients) of (|jtt'^||j^i^jjj, (o£:(it'')). Since 

(n'5^ y^(i|^j~(M^^y'+~a^^2(M^)(j converges almost everywhere to u\^dii {u\) + ai2(n2), and since S^{uP) converges to 
S{u^), the classical weak estimate ensures that 



|MiV<iii("i) +ai2('"2)| 
for some polynomial function P with positive coefficients. 



<P(|1uO||li,^„^(u»)) 



□ 
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